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Abstract 

A nonperturbative model for the QCD invariant charge, which contains no low-energy 
unphysical singularities and possesses an elevated higher loop corrections stability, is 
developed in the framework of potential approach. The static quark-antiquark poten- 
tial is constructed by making use of the proposed model for the strong running coupling. 
The obtained result coincides with the perturbative potential at small distances and 
agrees with relevant lattice simulation data in the nonperturbative physically-relevant 
region. The developed model yields a reasonable value of the QCD scale parameter, 
which is consistent with its previous estimations obtained within potential approach. 

PACS: 11.15.Tk; 12.38.Aw; 12.39.Pn 
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M ' 1 Introduction 

S. 

Theoretical description of hadron dynamics at large distances remains a crucial chal- 
lenge of elementary particle physics for a long time. The asymptotic freedom of Quan- 
tum Chromodynamics (QCD) allows one to apply perturbation theory to the study of 
some "short-range" processes, for example, the high-energy hadronic reactions. How- 
ever, the study of many phenomena related to the "long-range" dynamics (such as 
confinement of quarks, structure of the QCD vacuum, etc.) can be performed only 
within nonperturbative methods. 

In general, there is a variety of nonperturbative approaches to handle the strong 
interaction processes at low energies. For example, one can gain some hints about 
the hadron dynamics in infrared domain from lattice simulations pQ, string models [2], 
AdS/CFT methods [2111], sum rules [5JE], dispersive (or analytic) approach to QCD [7- 
E], bag models [ID] , potential models [TTHTB] . In what follows we shall employ the latter 
approach, that involves the construction of the QCD invariant charge which satisfies 
certain nonperturbative requirements. 
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The objective of this paper is to develop a model for the QCD invariant charge 
in the framework of potential approach. It is also of a primary interest to apply 
the proposed model to the construction of the static quark-ant iquark potential and 
compare the obtained result with relevant lattice simulation data. 

The layout of the paper is as follows. In Sect. |2] the model for the strong run- 
ning coupling is formulated and its properties are discussed. Section [3] contains the 
construction of the potential of quark-antiquark interaction and its comparison with 
lattice data. In Conclusions (Sect. H]) the obtained results are summarized and further 
studies within the approach on hand are outlined. Appendix |A] contains the explicit 
expressions for the perturbative QCD /3-function and strong running coupling up to 
the four-loop level. A brief description of the multi-valued Lambert IV-function is 
given in App. [B] 



2 A nonperturbative strong running coupling 

As it has been mentioned in the Introduction, the asymptotic freedom enables one to 
study the strong interaction processes at high energies within perturbative approach. 
However, the low-energy hadron dynamics entirely remains beyond the applicability 
range of perturbation theory. In what follows we shall adhere the so-called poten- 
tial approach [TTHT6] to the construction of the nonperturbative model for the QCD 
invariant charge. 

Specifically, in accordance with the basic idea of the potential approach, we shall 
construct the strong running coupling ct(Q 2 ) that coincides with perturbative QCD 
invariant charge in the ultraviolet domain 

a(Q 2 )~a peit (Q 2 ), Q 2 -> oo, (1) 

and meets the requirement of the infrared enhancement^ 

Here a pert (Q 2 ) is the perturbative strong running coupling (see App. IA"]) . f3 = 11 — 
2rif/3 denotes the first coefficient of the /3-function perturbative expansion, n { stands 
for the number of active quarks, Q 2 > is the spacelike momentum transferred, and 
A denotes the QCD scale parameter. In terms of the renormalization group /3-function 

d In a(a 2 ) _ . , 

-n V = fl (3) 

the conditions (JT]) and (T5]) can be equivalently rewritten as 

P»^/3 pert (a), a^0 + , (4) 

and 

15(a) ~ —1, a — Y oo, (5) 



^dt is worth noting here that the low-energy behavior of the strong running coupling ^ corresponds 
to the linearly rising at large distances static quark-antiquark potential in the framework of the one- 
gluon exchange model (see discussion of this issue in reviews |14H16j ). 
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respectively. In these equations a{Q 2 ) = a(Q 2 )(3o/(4:7i) stands for the so-called "cou- 
plant" . 

One of the possible expressions for the /3-function that satisfies conditions (pEJ) 
and (jSJ) read^l 



where pfl rt {a) is the £-loop perturbative QCD /3-function f ]A.2j) and B n = (3 n / '[3q +1 is 
the ratio of the /3-function expansion coefficients (see App. |A]). It is straightforward 
to verify that at any given loop level I the /3-function (j6j) satisfies the aforementioned 
condition of linear confinement (J5J): 



a V fl2 / 
where 

' 1, £= 1 



is the combination of the /3-function expansion coefficients. At the same time, for 
small values of the couplant a the /3-function fl6]) coincides with the corresponding 
perturbative result ( ]A.2j) up to the uncontrollable at £-loop level corrections: 

0®(a) ~ ^l(a) + 0(a m ), a 0+. (9) 

Additionally, /3-function ([6]) admits the explicit integration of the renormalization 
group equation (j3J) at the one-loop level and eventually leads to the strong running 
coupling a(Q 2 ) that contains no unphysical singularities at low energies and possesses 
an elevated stability with respect to the higher loop corrections. The plots of the func- 
tion (3^(a) dSJ) and the perturbative /3-function ( )A.2j) at £-loop level (£ = 1, ...,4) are 
presented in Figure [TJ 

For the beginning, let us consider the one-loop level (£ = 1). In this case the 
/3-function ([6]) takes a simple form, namely, 

0W(a) = --^-. (10) 
1 + a 

The corresponding renormalization group equation for the QCD invariant charge a(Q 2 ) 
rGcicls 

rfln[a (1 V)] a (1 V) , . 

rfln/i 2 l + aW(M 2 )' 

After splitting the variables and integrating in finite limits, Eq. ( ITT1) takes the following 
form: 

1 



In 



a (1) (Q 2 




'12) 



2 Similar modifications of the /3-function have been proposed in, e.g., Refs. [TTJ[T3] , whereas an 
equivalent modification of the running coupling itself was considered in, e.g., Ref. |12j . 
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Figure 1: The /3-function f3^\a) © (solid curves) and the perturbative /3-function 
/3 { Jl t (a) (O (dashed curves). The presented functions correspond to rif = 3 active 
quarks. The numerical labels indicate the loop level. 



where 



A 2 = Ql^a^(Ql)ex P 



Air 



0o aW(Qg) 



(13) 



denotes the square of the one-loop scale parameter and Q$ is the normalization point. 
In turn, Eq. (fT2l can be solved explicitly in terms of the multi-valued Lambert W- 
functionE), namely, 

(3 W {z) A 2 

see also Refs. JTSJQIJ]. It is worthwhile to note here that only the principle branch 
of the Lambert VF-function, Wo(x), is physically meaningful herein. The plots of the 
couplant corresponding to the invariant charge (Tl4"|) . a^\Q 2 ) = 1/Wq(z), and the one- 
loop perturbative couplant (1A.8I) . a^] rt {Q 2 ) = l/\nz, are presented in Figure [2] (solid 
and dashed curves labeled "(1)", respectively). 

The obtained strong running coupling (1T4"]) contains no unphysical singularities in 
the physical domain Q 2 > 0. By making use of the expansions (IB. 21) and flB.3|) one can 
show that ct«(Q 2 ) dHJ) possesses the required infrared enhancement, 



« (1) (Q 2 ) - 4 i 

Po 



- + 1 + 0(z) 



Q 2 -> o+, 



(15) 



and tends to the perturbative result (1A.8j) in the ultraviolet domain, 



« {1) (Q 2 ) 



4vr f 



In z 



O 



InQnz) 



Q 2 ->■ oo. 



(16) 



3 The definition of the Lambert ^-function and its properties are described in Ref. [IT] and briefly 
overviewed in App. [B] 
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Figure 2: The £-loop couplant a^(Q 2 ) corresponding to the /3-function (solid 
curves) and the perturbative couplant Op er t(Q 2 ) (dashed curves). The presented func- 
tions are computed for rif = 3 active quarks and normalized to the value a(Q 2 ) = 1/2 
at Q 2 = 5Qq. The numerical labels indicate the loop level. 



It is worth mentioning here that the low-energy behavior of the QCD invariant charge 
similar to that of Eq. (fl5|) has also been discussed in Refs. [2QH27] . 

Let us proceed now to the higher loop levels (£ > 1). Here, the renormaliza- 
tion group equation with the /3-function (jSJ) can only be integrated numerically. 
Nonetheless, the asymptotic behavior of the £-loop QCD invariant charge a^(Q 2 ) 
corresponding to the /3-function (jB]) can be found explicitly. Specifically, at any loop 
level a^{Q 2 ) possesses an enhancement in the infrared asymptotic, 



47T 



1 



+ C e + 0(z) 



Q 2 ^0 



;i7) 



and coincides with perturbative strong running coupling 0!^ Tt (Q 2 ) (see App. |Aj in the 
ultraviolet asymptotic, 



Q 2 ->■ oo. 



The constant Cg in Eq. ( TTT1) is defined in Eq. (jSJ). It is worth noting here that the last 
term of Eq. (fl8l) constitutes a correction, which is not controllable within perturbative 
approach at the £-loop level. The plots of the couplants a^(Q 2 ) and a p ert(Q 2 ) a ^ ^-loop 
level (£ = 1, 4) are presented in FigureO As one can infer from this figure, the QCD 
invariant charge a^(Q 2 ) possesses an elevated (with respect to perturbation theory) 
higher loop corrections stability in the intermediate energy range. It is worthwhile 
to note also that the proposed strong running coupling a(Q 2 ) is free of low-energy 
unphysical singularities. 
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3 Static quark— ant iquark potential 



Let us address now the construction of the static potential of quark-ant iquark inter- 
action V(r). In the framework of potential approach V(r) is related to the strong 
running coupling a(Q 2 ), which satisfies the aforementioned conditions ([[]) and (j2J), by 
the three-dimensional Fourier transformation 



V(r) 



16n f a(Q 2 ) exp(zQr) 



Q 2 



(2*) 



dQ, 



(19) 



see, for example, reviews [T4T[T5] and references therein. Strictly speaking, this defini- 
tion of the potential is justified for small distances (r < 0.1 fm) only. For instance, the 
lowest-lying bound states of heavy quarks can be described by employing the pertur- 
bativeB QCD [28j . However, at large distances (r > 0.5 fm), which play the crucial role 
in hadron spectroscopy, the perturbative approach becomes inapplicable due to the 
infrared unphysical singularities (such as the Landau pole) of the strong running cou- 
pling o> t {Q 2 ). Nonetheless, the potential (JT9~j) . being constructed with the invariant 
charge a(Q 2 ), which contains no unphysical singularities and satisfies requirements (JTJ) 
and has proved to be successful in the description of both heavy-quark and light- 
quark systems (see papers [TTHT3"] and reviews [TITHo] for the details). 

In this paper, for the construction of the static potential of quark-ant iquark inter- 
action we shall use the invariant charge (1141) . After integration over angular variables, 
Eq. (fl9l) in this case takes the form 



V(r) 



32 



3/?c 



oo 



(20) 



where p = Q/A, R = rA, and a(Q 2 ) = 1/W (Q 2 /A 2 ). The integral ([20]) diverges at the 
lower limit, that is a common feature of the models of such kind (a detailed discussion 
of this issue can be found in Sect. 7 and App. C of Ref. [15J). To regularize Eq. ( 1201) it is 
convenient to split the couplant f[T4"l) into singular and regular parts (see also Eq. ( TT5"|) ): 



a(Q 2 



where z = Q 2 / A 2 and 



ai(z) 



1 

z ' 



aAz) + a 2 u) 



a 2 (z) 



W n (z) 



1 

z 



Then, the potential ( |20l) takes the following form: 



57T 



V(r) = —A 

6p 



U 1 (R) + U 2 {R) 



where 



1 4 f | I" \S1M 





% = 1,2 



(21) 



(22) 



(23) 



(24) 



4 The leading short-distance nonperturbative effect due to the gluon condensate has also been 
accounted for in Ref. [251. 
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are the dimensionless parts of the potential and x = Qr. The function U\{R) diverges 
and requires regularization, whereas ^(-R) is regular and can be computed numerically. 

To regularize function U\(R) we shall employ the method similar to that of used 
in Ref. [20]. Specifically, in terms of auxiliary function 

oo 

I{t) = J X* sin x fix (25) 
o 

the singular part of the potential ( 1231) reads 

UJR) = I(-3)R. (26) 

7T 

The integral on the right hand side of Eq. ( 125]) exists for < |Re (t + 1)| < 1 (see, for 
example, Ref. [29]). Nonetheless, Eq. ( 1251) can be analytically continued to the entir^ 
complex t-plane. This continuation is given by 

J W = V ^ f + t {*\ (27) 

w r((i-t)/2) 

and plays the role of regularization of Eq. (|26|) . see also Refs. [U|20]. In this case 
I (-3) = -tt/4, that results in U^R) = R. 

Thus, the static quark-antiquark potential (12"U|) takes the following form: 



V(r) = V + -^A 




X 2 \ SIM 



R — / «2 d,X 



R 2 



x 



R = Ar, (21 



where Vo is an additive self-energy constant and a,2(z) is given by Eq. (1221 . At small 
distances this potential possesses the standard behavior determined by the asymptotic 
freedom 

^^jk^R^R' r ^°' (29) 



whereas at large distances potential (1281) proves to be linearly rising 

V(r)~^AR, r^oo, (30) 

implying the confinement of quarks. It is straightforward to verify that the poten- 
tial (128p satisfies also the so-called concavity condition 

dV(r) d 2 V(r) , , 

dr dr* 

which is a general property of the gauge theories (see Ref. [30] for the details). 

Figure |3] presents the quark-antiquark potential (|28|) . relevant lattice simulation 
data [31J, Cornell potential [32J, and Richardson's potential [12J. Equation ( 1281) has 
been fitted to the lattice data [21] by making use of the least square method, A and V 



5 Except for the singular points of the right hand side of Eq. (|27|) . such as t = —2N, with N 
1, 2, 3, ... being a natural number. 
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Figure 3: The quark-antiquark potential V(r) (Eq. fl28|) . solid curve) and the lattice 
simulation data (Ref. [31], symbols "O"). The Cornell potential [32] is denoted by sym- 
bols whereas Richardson's potential [12] is shown by symbols "A". The dashed 
curve denotes the relevant perturbative result (1291) . The values of the parameters: 
A = 375 MeV, V = 315 MeV, n f = 3. 

being the varied parameters. The estimation of the scale parameter in the course of this 
comparison gives A = (375 ±40) MeV (this value corresponds to the one-loop level with 
nf = 3 active quarks). The obtained value of A agrees with previous estimations [T1T[T3] 
of the QCD scale parameter within potential approach. As one can infer from Figure [3l 
in the region r < 0.05 fm the derived potential (T28]) coincides with the perturbativej 
result (I29p . At the same time, V(r) ( 12 8 p is in a good agreement with both Cornell [32] 
and Richardson's [T2] potentials. Besides, in the nonperturbative physically-relevant 
range 0.3 fm < r < 1.2 fm, in which the average quark separations \J (r 2 ) for quarkonia 
sits [34], the obtained potential ( 128|) reproduces the lattice data [31] fairly well. 

4 Conclusions 

A nonperturbative model for the QCD invariant charge is developed in the framework 
of potential approach. The proposed strong running coupling is free of low-energy un- 
physical singularities and embodies the asymptotic freedom with infrared enhancement 
in a single expression. The model on hand possesses an elevated (with respect to per- 
turbation theory) higher loop corrections stability in the intermediate energy range. 
The static quark-antiquark potential is constructed by making use of the proposed 
model for the strong running coupling. The obtained result coincides with the per- 
turbative potential at small distances and agrees with relevant lattice simulation data 
in the nonperturbative physically-relevant region. The developed model yields a rea- 
sonable value of the QCD scale parameter, which agrees with its previous estimations 

6 Thc applicability range of the perturbative quark-antiquark potential has been discussed in 

Refs. rrji^]. 
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obtained in the framework of potential approach. 

In further studies it would undoubtedly be interesting to compute the meson spec- 
trum by making use of the derived quark-ant iquark potential as well as to apply the 
developed model for the QCD invariant charge to the analysis of other strong interac- 
tion processes. 
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A The QCD function and running coupling 
within perturbation theory 

The QCD invariant charge a(Q 2 ) is the solution of the renormalization group equation 



d In [a(/i 2 )] 
d In u 2 



/9(a), 



(A.l) 



where a{Q 2 ) = a(Q 2 )f3 / (4ti) is the so-called "couplant". For small values of the 
running coupling the right hand side of Eq. 1 1A.1I) is usually approximated by the 
power series, namely, 



e-i 

/J(a)~/3&(a) = -£i^V 



n=0 



n+1 



(A.2) 



where I = 1, 2, 3, ... denotes the loop level and B n = f3 n / (3q +1 is the ratio of the QCD 
/3-function perturbative expansion coefficients: 



Pi 



11 Tlf, 

3 ' 
38 

102 n { , 

3 

2857 5033 



2 1 
149753 



7lf + 



325 



54 



■ n 



6 



+ 3564C(3) 



f ) 

1078361 + 6508 ' 



162 



27 



+ 



50065 6472 > , ,' 



2 , !093 ^3 



n f + 



729 



(A.3) 
(A.4) 
(A.5) 



(A.6) 



In these equations nf stands for the number of active quarks and ((x) denotes the 
Riemann ^-function, £(3) ~ 1.202 (see, for example, Ref. [29]). The one- and two- 
loop coefficients (/3 and {3\) are scheme-independent, whereas the expressions given 
for 02 and /3 3 are calculated in the MS subtraction scheme (see papers [35H38] and 
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references therein for the details). The plots of the /3-function ( 1A.2j) are presented in 
Figured] (dashed curves). 

The renormalization group equation corresponding to the perturbative /3-function 

d\n[af ert {^)] _ (fl 
can be solved explicitly at one- and two-loop levels (£ — 1,2), namely, 



a pert(Q ) 



1 



In z 



(A.8) 



z (2) 

"pert 



(Q 2 



l + W_i{-exp -(l + B^lnz) } 



In 2; 



ln(ln z) 
ln 2 z 



Q 2 



00, 



(A.9) 



see also Ref. j3H]. In Eq. (1A.9j) M4(x) stands for the so-called Lambert W^-function 
(see App. |B]for the details). Starting from the three-loop level, the exact solution to 
the perturbative renormalization group equation (IA.7j) can not be expressed in terms 
of known functions. Nonetheless, for I > 3 Eq. (1A.7|) can be solved iteratively, that 
eventually leads to 



a pert{Q' ) 



1 



In ,2 



ln(lnz) 
In 2 z 



ln d z 



B\ ln 2 (lnz) -ln(lnz) - 1 



B, 



(A.IO) 



a pert{Q 



Inz 



Bi 



\n(lnz) 
In 2 z 



ln 3 z 



B\ ln 2 (lnz) -ln(ln^) - 1 



Br 



ln 4 z 



B\ 



ln 3 (lnz) + - In 2 (In z) 



+21n(lnz) 



3BiB 2 \n(\nz) + -B 3 



(A.ll) 



It is worth noting here that the presented in Figure [2] plots of the £-loop perturbative 
couplant dp e rt(Q 2 ) (dashed curves) correspond to exact solutions of the renormalization 
group equation (IA.7j) . 
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B The Lambert W— function 



The Lambert J^-function is defined as the multi-valued function Wk(x), which satisfies 
the equation 



W k (x) exp 



WJx) 



(B.l) 



where k denotes the branch index. Only two branches of this function, namely, Wq(x) 
(the principal branch) and W-i(x), take real values (see Figure Hj), whereas the other 
branches take imaginary values. For the branches Wq(x) and W-\(x) of the Lambert 
W-function the following expansions hold: 



W (x) 

Wo(x) 
W_i(-z) 

h x 



Inx — ln(lnx) + O 



ln(lnx) 



lnx 

x-x 2 + 0(x 3 ), x -» 0, 
lnx + 0(\i\ | lnx|), x — > + , 

-1- V2ex + 0(x), x-^0- 



x — > oo, 



(B.2) 

(B.3) 
(B.4) 

(B.5) 



where e ~ 2.71828 denotes the base of natural logarithm. The detailed description of 
the Lambert ^-function and its properties can be found in Ref. [17] . 




Figure 4: Two real branches of the Lambert jy-function fIB.lj) : Wq{x) (solid curve) 
and W-i(x) (dashed curve). 
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